Kinetic Roughening of Ion-Sputtered Pd(OOl) Surface: Beyond the Kuramoto-Sivashinsky Model 
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We investigate the kinetic rougliening of Ar+ ion-sputtered Pd(OOl) surface both experimentally and the- 
oretically. In situ real-time x-ray reflectivity and in situ scanning tunneling microscopy show that nanoscale 
adatom islands form and grow with increasing sputter time t. Surface roughness, W{t), and lateral correlation 
length, t,{t), follows the scaling laws, W{t) ~ fP and t,{t) ~ with the exponents (3 ~ 0.20 and 1/z ~ 0.20, 
for ion beam energy £ = 0.5 keV, which is inconsistent with the prediction of the Kuramoto-Sivashinsky (KS) 
model. We thereby extend the KS model by applying the Sigmund theory of sputter erosion to the higher order, 
o(V^,/!^), where h is surface height, and derive a new term of the form V^{Vh)^ which plays an indispensable 
role in describing the observed morphological evolution of the sputtered surface. 

PACS numbers: 68.55.-a, 05.45.-a, 64.60.Cn, 79.20.Rf 



Recently, the observation of ordered nanostructures such 
as ripples and two-dimensional patterns on ion-sputtered sur- 
faces has attracted much attention due to the demonstration of 
the possibility of fabrication of ordered nanoscale^ structures 
in a relatively easy and affordable way 
Such experimental results have motivated extensive theoreti- 
cal investigations aiming to understand the mechanism of the 
morphological evolution of ion-sputtered surfaces. A linear 
model, proposed by Bradley and Harper (BH) |8], has been 
successful in predicting the formation of the ripple structure. 
The wavelength, orientation and amplitude of the ripples can 
be predicted in terms of experimental parameters such as the 
incident angle of the ion beam and substrate temperature |9]. 
The BH theory, however, fails to explain a number of ex- 
perimental observations such as the saturation of the ripple 
amplitude [HI [T3l . or the appearance of kinetic rough- 
ening Il4ll . To remedy such shortcomings, the noisy 
Kuramoto-Sivashinsky (KS) equation 1 15, 16] was introduced 
based on the Sigmund theory of sputter erosion 1 17]. In addi- 
tion to the linear terms of the BH model, it contains a nonlin- 
ear term proportional to (V/i)^, known as the Kardar-Parisi- 
Zhang (KPZ) term ]18], where h is surface height. Due to 
the nonlinear term, the surface roughness (or the ripple ampli- 
tude), which was growing exponentially with increasing sput- 
ter time in the linear model, changes to the type following a 
power law and eventually saturates to a constant value ] 19]. 
Although the KS model seems to be successful in offering 
insights for understanding the nonlinear behavior of sputter- 
eroded surfaces, it has not been convincing yet because the 
detailed properties of the kinetic roughening predicted by the 
KS equation have not been fully tested experimentally. 

Kinetic roughening behavior is described by the scal- 
ing theory L20.I . The surface roughness of the sam- 



ple, L X L in size, at sputter time t is defined as W = 
\/{^/L^)Lr[h{r,t)-hf, where h = (l/L^) ^^/,(r,f). It 
follows the scaling relation, W{L,t) ^ L°'g{t / L°'l^), where 
g(u) ~ for M ^ 1 and g[u ^ oo) ^ constant. The rough- 
ness and growth exponents a and p are related via scaling 
relations to give the dynamic exponent, z = a/p, which deter- 
mines the scaling of the saturation time with the system size L. 
Below, we will deal with the inverse dynamic exponent, 1 /z, 
called the coarsening exponent, in analyzing the experimental 
results. 

In this Letter, we study the kinetic roughening of sputter- 
eroded Pd(OOl) surface via novel experimental techniques 
combined with the stochastic continuum theory. Both in situ 
real-time X-ray reflectivity (XRR) and in situ scanning tunnel- 
ing microscopy (STM) show that nanoscale islands are formed 
on surface, evolving with increasing sputter time t following 
the scaling function W(L,f). In particular, we obtain p ~ 0.2 
and \/z — 0.2 for ion beam energy e = 0.5 keV, which are 
not in agreement with the values predicted by the KS model. 
To resolve this inconsistency, we investigate theoretically the 
erosion process by applying the Sigmund theory to higher or- 
der, finding that the relevant higher order term is of the form 
V^(V/!)^, referred to as the conserved KPZ (cKPZ) term. Us- 
ing this new nonlinear term, we can explain the kinetic rough- 
ening on the sputtered Pd(OOl) surface successfully. 

Experiment. — A Pd(OOl) sample was cleaned by several 
cycles of Ar+ ion sputtering around 300 K and annealing up 
to 920 K. The clean Pd(OOl) surface exhibited an average ter- 
race size of about 3000 A and little surface modulation or de- 
fects by contaminants as judged by STM images. The initial 
surface roughness of the sample was less than 2 A and the 
miscut from surface normal direction was 0.3°, as determined 
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FIG. 1: Specular x-ray reflectivity (symbols) as a function of the 
momentum transfer, q,, for increasing sputter times with e = 0.5 keV 
and / = 2.0 x lO'"' ions/cm^/sec. Sputtering is made normal to a 
sample at room temperature. The solid lines are the theoretical ones 
according to the Parratt formula. 
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FIG. 2: The surface roughness IV as a function of sputter time / for 
different ion energies, e = 0.5(O)> l-5(n), and 2.0(A) keV, but with 
a fixed ionic flux, / = 2.0 x lO'-' ions/cm^/sec. Inset shows that the 
roughness evolution for a given e follows the same fluence depen- 
dence, irrespective of /. 



respectively by X-ray reflectivity (XRR) and X-ray diffraction 
spot-profile analysis. Sputtering experiments were performed 
for Pd(OOl) around 300 K with Ar+ ion beam incident nor- 
mally to the sample surface. To avoid possible contamination 
during sputtering, fresh Ar gas (purity of 99.999%) was con- 
tinuously made to flow through the chamber with Ar partial 
pressure maintained around 2.0 x 10^^ Torr 

The morphological evolution of the Pd(OOl) during sputter- 
ing was observed in situ by both XRR and STM. The XRR 
experiment was performed in a custom-designed UHV cham- 
ber at the 5C2 K-JIST beamline of the Pohang Light Source 
in Korea. The base pressure of the chamber was kept below 
mid- 10^^" Torr. The incident X-rays were vertically focused 
by a mirror to 1 mm in size, and a double-bounce Si(lll) 
monochromator was employed to both focus the beam in the 
horizontal direction and monochromatize the incident X-ray 
to a wavelength of 1.24 A which corresponded to the X-ray 




FIG. 3: X-ray diffraction intensity along [110] for different q^. The 
positions of satellite peaks shift to the larger g|| for the larger q,, with 
a fixed angle to the surface normal, indicating that the spots originate 
from a facet formed by sputtering. Inset: Plot of the positions of 
satellite peaks, , as a function of q,. From the angle the facet peaks 
form with respect to the surface normal, we identify the facet planes 
to be (117) for e = 0.5 keV and (115) for e = 2.0 keV. 



energy of 10 keV. The X-ray reflectivity was measured in real 
time while Pd(OOl) was sputtered. Incident Ar+ beam en- 
ergy, 8, ranged from 0.5 to 2.0 keV, and ion beam flux, /, 
measured from the ion current collected at the sample, ranged 
from 0.5 x lO'^ to 2.5 x lO'-' ions/cm^/sec. Sputtering started 
with the sample at room temperature measured by a thermo- 
couple directly attached to the back-end of the sample. Af- 
ter sputtering for more than 10 hours, the sample temperature 
rose only less than 30 K. 

To observe the sputter-induced morphological evolution in 
real-space, a commercial STM (Omicron VT-STM) was also 
used. The base pressure was kept in the low 10^^' Torr. STM 
tips were electro-chemically etched poly crystalline W- wires 
that were annealed and sputtered in a UHV chamber In all 
the STM measurements, the tip-sample voltage was kept as 
80 mV, and the tunneling current was at 1 nA with the sample 
at room temperature. Sputtering was performed with 8 = 0.5 
keV and / = 0.5 x 10^^ ions/cm^/sec. During sputtering, the 
sample temperature was around 300 K. Sputtering and STM 
measurement were alternated, and the sputter time referred 
to in the present work is the total sputter time determined by 
summing all the previous sputtering periods. X-ray reflectiv- 
ity indicated that there was no significant change in the sur- 
face roughness for 10 hours after the interruption of sputter- 
ing. Thus, the inteiTuption of sputtering for the STM measure- 
ment was not expected to result in additional relaxation of the 
sputtered surface. 

Experimental results. — In Born approximation which holds 
above critical angle, the specular reflectivity decreases as 
exp(— W^<7^) with increasing |21]. Fig. 1 shows that the 
reflectivity decays more steeply as the sputtering proceeds, 
which indicates that W{t) increases with t. Quantitatively 
W{t) is determined by fitting the experimental data accord- 
ing to the Parratt formalism L21i. i22.] . Fig. 2 shows W{t) thus 
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FIG. 4: STM images for various sputtering times whicii are 200 x 
200 nm^ in size. The insets at right bottom of each image show the 
height-height correlation function obtained from the corresponding 
surface profile. 
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FIG. 5: Plot of the lateral correlation length ^(O) ''id the surface 
roughness lV(n) obtained from the STM images as a function of 
sputter time t. 
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obtained, as a function of t for several £ but with a fixed /. The 
linear behavior in double-logarithmic scale implies power-law 
behavior, IV (f) ~ fl^. We obtained that p ~ 0.20 ±0.02 for e = 
0.5 keV, p ~ 0.23 ± 0.02 for e = 1 .5 keV, and P ~ 0.25 ± 0.02 
for e = 2.0 keV. 

The inset of Fig. 2 summarizes the surface roughness, W{t) 
as a function of fluence (/ x f) for several different /s, but 
with a fixed e. We find that the roughness evolves showing 
the same fluence dependence, irrespective of /. Such behav- 
ior implies that the surface diffusion is not induced by ther- 
mal activation, but rather by the surface sputtering under the 
present experimental condition ll23il . 

We also study the facet formation on the sputtered surface 
by tracing the extra-diffraction peaks or satellites around the 
diffraction peak corresponding to the sample surface, now 
(001) plane, because each facet plane defines diffraction rods 
normal to it. In the inset of Fig. 3, the diffraction rods defined 
by the satellite peaks make well-defined angles with respect 
to that of (001) plane as much as the facet plane does to the 
(001) plane . From the observed angles, we identify the facet 
formed on sputtering by 0.5 keV ion beam to be (117) and that 
by2.0keV beam to be (115). 

The morphological evolution of the sputtered Pd(OOl) sur- 
face were also investigated by an STM. Fig. 4 shows the STM 
images of the sputtered surface at different (total) sputtering 
time under the same experimental condition, e — 0.5 keV and 
/ = 0.5 X 10^^ ions/cm^. The shape of the island edge looks 
somewhat irregular which may be due to random shot noise 
or the sputter-induced diffusion. This reaffirms that thermal 
smoothening is inefficient in the present experimental condi- 
tions as suggested by the flux-independent scaling behavior. 
The insets of Fig. 4 show the height-height correlation func- 
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FIG. 6: Comparison of the KPZ and the cKPZ terms for a prototypi- 
cal profile of eroded surface. 



tion, G(r,f) = {{h{r,t) -/i(0,0))^), that reveals lateral order 
developing on the surface as sputtering proceeds. We can see 
that the island size grow with increasing sputter time. From 
the height-height correlation function, we determine the lat- 
eral correlation length, ^, as the position giving the first-order 
maximum 1I24I1 . and plot it as a function of t in Fig. 5. It 
is found that ^ f^/^ with the exponent l/z ~ 0.20 ±0.03. 
The growth exponent determined by the STM measurement is 
P ~ 0.20 ±0.03, in agreement with that of the XRR measure- 
ment. 

Theory. — To elucidate the kinetic roughening, we system- 
atically expand the Sigmund theory to higher order than that 
where the KS equation is derived, and obtain the dynamic 
equation, 

- ^ = - 1 - V V^/; - DV^h + 'ki(Vhf + A.2 V2 (Vhf +\\{x,y,t) , 
c dt 

(1) 

where c is the surface recession rate given by the material 
constant times the ionic flux, v the effective surface tension 
generated by the erosion process, D the ion-induced effective 
diffusion constant |23|, Xi and A-2 the tilt-dependent erosion 
rate, and r| an uncorrected white noise with zero mean, mim- 
icking the randomness resulting from the stochastic nature of 
ion arrival at the surface. The coefficients v, D, and Xi in 
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Eq. Q were derived in terms of experimental parameters such 
as the ion beam energy e, the flux /, the incident angle 9, the 
penetration depth a, and the cascading sizes in transverse and 
longitudinal directions, /j and o respectively 1 16|. The newly 
derived cKPZ coefficient is, under normal incidence, given by 
A.2 = A'V2 + (3/8)(a(/o)'*(o^ -fl^). Using the trim Monte 
Carlo algorithm |25|, all the coefficients in Eq. Q are de- 
termined numerically for our experimental parameters, which 
are tabulated in Table I. In Fig.|6l is shown the estimation of 
the KPZ and cKPZ terms, for a typical profile of the eroded 
surface. We find that the cKPZ term should be more relevant 
to the morphological evolution than the KPZ term in a finite 
system. The contributions from other terms allowed by sym- 
metry and of order o(V'*,/i^) such as (V^/z)^, are found to be 
negligible. Then, for e — 0.5 keV, the extended KS equation is 
reduced to the cKPZ equation that gives the growth exponent 
P ~ 0.20 1 26| for two dimensions, in agreement with the mea- 
sured value for e — 0.5 keV. As the ion energy increases, the 
coefficient D increases very rapidly in comparison with the 
magnitudes of the A-i and A,2 as seen in the Table I. Thus, for 
the e = 2.0 keV, the erosion process is mainly governed by the 
so-called Mullins term, DV^/z fxf]. In this case, it is known 
that the growth exponent p ~ 0.25 for two dimensions, again 
consistent with the experimental value for e = 2.0 keV. The 
sputter-induced palladium atoms form adatom islands with 
well-defined facets as observed in Fig. 3. Then, a = 1 and 
the coarsening exponent, 1 /z, should be equal to P following 
the scaling relation 1/z = p/a. This prediction is confirmed 
in the present experiment, by the observation of 1/z = p = 0.2 
for e = 0.5 keV. 
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TABLE I: Numerical estimations of the coefficients in Eq.Q under 
the experimental conditions we performed. 

Conclusions. — We have studied the kinetic roughening of 
the sputter-eroded Pd(OOl) surface both experimentally and 
theoretically. The experimental data suggest that the KS con- 
tinuum model is not relevant to the morphological evolution 
of the Pd surface. Instead, the KS model extended by the in- 
clusion of the cKPZ term we derived from the Sigmund theory 
could properly describe the kinetic roughening behavior of the 
sputtered Pd(OOl) surface. The model, at the same time, may 
shed light on recent experimental realization of self-organized 
surface nanostructures induced by ion sputtering. For exam- 
ple, the parameter space for nanoscale island formation pre- 
dicted by the KS model, which seems too narrow to be real- 
istic as compared with experimental observations [28,1 . can be 
enhanced. 
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